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an effective potential V e f. Using a finite differences method and the Crank-Nicholson scheme, in a code 
implemented in the program OCTAVE, we solve the corresponding time dependent Schrodinger equation 
and obtain the time evolution of a initial wave packet. This wave packet satisfies appropriate boundary 
conditions. The calculation of the tunneling probabilities shows that the universe may emerge from the 
Planck era to an inflationary phase. It also shows that, the tunneling probability is a function of the mean 
energy of the initial wave packet and of two parameters related to the Chaplygin gas. We also show a 
comparison between these results and those obtained by the WKB approximation. 
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I. INTRODUCTION 

One of the main purposes of quantum cosmology is to fix the initial conditions that will determine the later behavior 
of the Universe In quantum cosmology, the Universe is treated as a quantum system. Since the cosmological 
scenario is based on general relativity, the dynamical variables are directly related to the geometry of the space-time. 
This implies a functional equation, the Wheeler-DeWitt equation, on all possible geometries [2j. Such functional 
equation has in general no exact solution. A possible way to circumvent this technical limitation is to freeze out an 
infinite number of degrees of freedom, remaining at the end with what is called a minisuperspace, a configuration 
space with few degrees of freedom, which may admit exact solutions. Even if such procedure may be seen as a drastic 
limitation of the original problem, it allows at least to extract some precise predictions concerning the evolution of the 
early Universe. One of the hopes in such a program is to obtain initial conditions for the classical evolution leading 
to an inflationary expansion of the Universe. 

Yet, the limitation to a minisuperspace approach does not eliminate other difficulties in quantum cosmology. The 
Wheeler-DeWitt equation has no explicit time variable [3( . This is a consequence of the fact that the Einstein- Hilbert 
action represents a constrained system which is invariant by time reparametrization. Moreover, the limitation to 
a few number of degrees of freedom, the restriction to the minisuperspace, does not assure that the final equation 
will admit closed solutions. Frequently, it is necessary to look for approximative solutions [1] using, for example, 
perturbative methods. The usual WKB method of ordinary quantum mechanics is one of the possible methods to 
extract predictions of a quantum cosmological system [5| . It is a quite attractive procedure since it allows to introduce 
a time variable which is, otherwise, absent. 

In some specific situations it is possible to recover the notion of a time variable in the minisuperspace approach. 
One example is when gravity is coupled to a perfect fluid, which is expressed through the Schutz's variable [q]: the 
conjugate momentum associated with the fluid's variable appears linearly in the Lagrangian, which allows to recover 
a genuine Schrodinger equation Since now the quantum system has an explicit dynamics, it is possible to study 
the emergence of the classical Universe in a primordial era. An important question can be addressed in this way: 
what are the conditions to have an initial inflationary phase in the evolution of the Universe? 

We will study this problem in this paper. Gravity will be coupled to a radiative perfect fluid, whose dynamic 
variables will be connected with the time evolution of the quantum system. Aside the radiative fluid, the Chaplygin 
gas will be included in this scenario. The Chaplygin gas is an exotic fluid exhibiting negative pressure which depends 
on the inverse of the density |8(. The interest on the Chaplygin gas fluid has been, until now, mainly connected with the 
existence of a dark energy component in the Universe today Q. In this sense, it leads to a scenario competitive with 
the AC DM model. However, it can be of interest for the early Universe, and this due to a main reason: the Chaplygin 
gas can be obtained from the string Nambu-Goto action re-expressed in the light cone coordinate [1 01 ] - Since string 
theory manifests, in principle, its main features in very high energy levels, as it occurs in the early Universe, exactly 
when the quantum effects must be relevant, a quantum cosmological study of the Chaplygin gas is of course relevant. 
Moreover, the Chaplygin gas admits a de Sitter asymptotically phase, that is, it tends to a cosmological constant, 
being in one sense a much richer structure. In that sense, one may consider the present paper a generalization of a 
previous paper made by some of us, where instead of the Chaplygin gas we used a positive cosmological constant fllj | . 
For all these reasons, a quantum cosmological scenario using radiative fluid and the Chaplygin gas is a very attractive 
configuration by its own. 

A quantum cosmological model with two fluids, as it will be treated here, implies considerable technical difficulties. 
In particular, no analytical solution can be expected. Hence, we are obliged to consider a perturbative analysis, like 
WKB [5], or a numerical integration of the equations. We will concentrate, here, in the latter possibility. It will be 
used a finite differences method and the Crank-Nicholson scheme, in a code implemented in the program OCTAVE. 
Using such program, we will compute the tunneling probability (TP) from ability (TP) from a quantum phase to 
the classical regime, the potential barrier being represented by the density of the Chaplygin gas in terms of the scale 
factor. When the system emerges from the potential barrier, it finds itself in an inflationary phase. The main results 
of our paper, besides the solution of the present quantum cosmological model, is related to the determination of the 
initial conditions for the classical evolution of the universe. Here, for the present model, we would gain information on 
what is the most probable amount of radiation in the initial evolution of the classical universe and the most probable 
values of the parameters of the Chaplygin Gas. If we take in account that one of those parameters tends to the 
cosmological constant (A), for great values of the scale factor, we would gain information on the most probable value 
of A. Another important result, derived here, is the fixation of the initial value of the scale factor and its conjugated 
momentum leading to a classical inflationary evolution of the scale factor. 

In fact, the quantum cosmology of a closed FRW model coupled to a Chaplygin gas have already been studied 
in the literature [T^. We may mention several differences between our work and Reference [T^. In their work, 
the authors considered the generalized Chaplygin gas as the only source of matter. Here, we consider the standard 
Chaplygin gas and radiation as our sources of matter. We use the Schutz formalism in order to introduce a time 
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variable and transform the Wheeler-DeWitt equation in a Schrodinger equation. Most importantly, we solve the 
Wheeler-DeWitt equation exactly using a numerical procedure, what was even suggested by the authors of Reference 
fl^ |. On the other hand, the authors of [13] use an approximated potential to describe the Chaplygin gas and solve 
the resulting Wheeler-DeWitt equation, to that approximated potential, in the WKB approximation. In despite of 
all those differences between the models and the treatments used here and in [12| , it is important to mention that we 
have obtained few qualitative agreements concerning the tunneling probability as a function of the parameters A and 
B. Their approximated, analytical expression for TP as a function of A, or the cosmological constant, grows with A 
for the boundary conditions of the tunneling wave-function [131 ]. Similarly, their approximated, analytical expression 
for TP as a function of B, grows with B for the boundary conditions of the tunneling wave- function [l3j]. Those 
results are in agreement with our numerical curves for TP as a function of A and as a function of B. 

This paper is organized as follows. In next section, the classical model is presented. In section 3, the classical model 
is quantized. In section 4, the wave packet describing the quantum system is constructed. In section 4 and 5, the 
tunneling effect is investigated and the initial conditions for inflation are settled out. The conclusions are presented 
in section 6. 



II. THE CLASSICAL DYNAMICS 



The homogeneous and isotropic Friedmann-Robertson- Walker (FRW) models with positive curvature of the spatial 
sections, radiation and Chaplygin gas, may be represented by a Hamiltonian with two degrees of freedom in the form 



H=—p 2 a + V ef (a)-p T , 



(1) 



where p a and px are, respectively, the moments canonically conjugated to the scale factor and to the variable that 
describes the perfect fluid, and V e f(a) is an effective potential, Figure [TJ. We are working in the unit system where 
G = 1 , c = 1 and Ti = 1 . V e f (a) is given by, 



y e/ ( a ) = 3a 2 -^-^+^, 



(2) 



where A and B are parameters associated to the Chaplygin gas. The Hamilton equations, 
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rule the classical dynamics of these models and are completely equivalent to the Einstein equations. The dot stands 
for ordinary differentiation relative to the conformal time r. Combining the first two equations in ([3]) we have, 
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(4) 



According to ((H) , the universe has an accelerated phase (d > 0) for big values of the scale factor (a 6 » B /A) and a 
de-accelerated one (a < 0) when a 6 << B/A. The solutions of (|4]) for the latter phase are described by trigonometric 
functions. 
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FIG. 1: Behavior of the effective potential V ef (a) for A = 0.001 and B = 0.001. 



III. THE QUANTIZATION OF THE MODEL 

The Planck era will be described by the quantization in the minisuperspace of that model, following the prescription 
proposed by Wheeler and De Witt [2[. This description is based on the wave function of the universe, which depends 
on the canonical variables a e T, 

* = ^(a,f) . (5) 

The quantization process begins by promoting the canonical moments p a and pr to operators and imposing that 
the Hamiltonian annihilates the wave function, 

Q Q 

Pa^-ifa > PT^-i-gf , H*(o,T) = 0. (6) 

Then, the quantum dynamics turns out to be ruled by the Wheeler-DeWitt equation, which in this case, assumes 
the form of a time dependent Schrodinger equation, 

where we impose the reparametrization r = —T. 

The operator H is self-adjoint [T3| in relation to the internal product, 



(*,*)= / daV(a,Ty$(a,T), (8) 
Jo 

if the wave functions are restricted to the set of those satisfying either \t(0, r) = or \I/'(0,t) = 0, where the prime 
/ means the partial derivative with respect to a. Here, we consider wave functions satisfying the former type of 
boundary condition and we also demand that they vanish when a goes to oo. 

The Wheeler-DeWitt equation (|7|) is a Schrodinger equation for a potential with a barrier. We solve it numerically 
using a finite difference procedure based on the Crank-Nicholson method [l6j . [l7j and implemented in the program 
GNU-OCTAVE [l8|]. Following the discussion in Ref. [llj . we also use, here, the norm conservation as the criterion 
to evaluate the reliability of our numerical calculations of the wave-packet's time evolution. Therefore, we have 
numerically calculated the norm of the wave packet for different times. The results thus obtained show that the norm 
is preserved. Numerically, one can only treat the tunneling from something process, where one gives a initial wave 
function with a given mean energy, very concentrated in a region next to a = 0. That initial condition fixes an energy 
for the radiation and the initial region where a may take values. Our choice for the initial wave function will be 
described in the following section. 
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IV. WAVE PACKETS 

As the initial condition, we have chosen the following wave function, 

where E m represents the mean kinetic energy of the wave packet associated to the energy of the radiation fluid. 
Besides, this initial condition is normalized as follows J °° \^(a,0)\ 2 da = 1. 

The portion of the wave function that tunnels the potential barrier, propagates to infinity in the positive scale 
factor direction, almost like a free particle, as time goes to infinity. However, we must specify a limit, in the scale 
factor direction, in order to perform the numerical integration of the Schrodinger equation. In this case, this value 
was fixed in a = 30. The behavior of these wave packets and their time evolution show that they are well defined in 
the whole space, even when a goes to zero. 

As an example, we fix the values of the system parameters A = 0.001 and B = 0.001. We choose the initial wave 
packet with mean energy E m = 220, which gives an initial amplitude of 5.901650048. This wave packet will evolve in 
time, obeying the Schrodinger equation (J7J), and reaches the potential barrier at the point a\ = 11.41476507. In this 
case, the potential barrier has a maximum value V ma x — 223.5282023 at a max — 12.20732711. The figures show the 
time evolution of these wave packets ( |W(a, 0)| 2 ), which indicates the possibility of the universe to have a posterior 
inflationary phase via a tunneling mechanism. With this choice for the parameters, the universe emerges at the right 
side of the potential barrier with the size — 12.95147904. 



0.5- 



0.4- 




a 



FIG. 2: \^(a,t m ax)\ 2 = p, for A = 0.001, B = 0.001 and E m — 220 at the moment t ma x when ^ reaches infinity, located at 
a = 30. 



V. THE TUNNELING EFFECT ON THE EARLY UNIVERSE AND THE INITIAL CONDITIONS FOR 

INFLATION 



A. Tunneling probability as a function of E m and comparison with WKB 

In order to evaluate the tunneling probability as a function of the mean energy E mi one needs to fix the values of 
A and B and repeat the calculations for different values of E m in *f?(a, 0). For each value of E m , the turning point 
associated to the potential barrier, to the left (ai) and to the right (0,2) are calculated. 

The probability to find the universe to the right side of the barrier is defined as [TTJ] , 
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TP 



|*(o, tmax)\ 2 da 
Jo°° \^( a > tmax)\ 2 da 



(10) 



where, as we have mentioned above, numerically oo has to be fixed to a maximum value of a. Here, we are working 
with a max = 30. It is useful to remember that the denominator of the equation (|10p is equal to 1, since the wave 
function is normalized to unity. 

In order to evaluate the tunneling probability as a function of the mean energy E m , one needs to fix the values of 
A and B and repeat the calculations for different values of E m in ^(a, 0). For each value of E m , the turning point 
associated to the potential barrier, to the left (a\) and to the right (02) are calculated. For all cases, we consider 
the situation where E m is smaller than the maximum value of the potential barrier. From that numerical study we 
conclude that the tunneling probability grows with E m for fixed A and B. As an example, we consider 60 values of 
the radiation energy for fixed A — 0.001 and B = 0.001. For this choice of A and B the potential barrier has its 
maximum value equal to 196.0471982. In order to study the tunneling process, we fix the mean energies of the various 
^(OjOj's eq. © to be smaller than that value. In table [HI in the appendix, we can see, among other quantities, the 
different values of the energy E m , TP, a\ and 02 for each energy. In figure [3j we see the tunneling probability as 
functions of E m , for this particular example. Due to the small values of some TP's, we plot the logarithms of the 
TP's against E m . 

Based on the detailed discussion made in Ref. [ill ], it is not difficult to understand under what circumstances the 
tunneling probability computed with the WKB wave- function (TPwkb) agrees with TP Eq. fTO]) . In order to see it, 
let us consider that, initially, we have an incident wave (\P 1) that reaches the potential barrier at 01. Then, part of 
is transmitted to 00 (vE't) and part is reflected (^ r)- In the present problem, we have an infinity potential wall at 
a = because the scale factor cannot be smaller than zero. It means that ^> r cannot go to —00, as was assumed in 
order to compute the TPwkb- Instead, ^r will reach the infinity potential wall at a = and will be entirely reflected 
back toward the potential barrier, giving rise to a new incident wave ((^ r)i). The new incident wave (^r)i reaches 
the potential barrier at a\ and is divided in two components. A reflected component which moves toward the infinity 
potential wall at a = (((^r)i)r) and a transmitted component which moves toward 00 (((\I> r)i)t)- ((^r)i)t will 
contribute a new amount to the already existing TP due to (^t)- Therefore, the only way it makes sense comparing 
TP with TPwkb is when we let the system evolve for a period of time (At) during which ^ r cannot be reflected 
at a = and come back to reach the potential barrier. It is clear by the shape of our potential that the greater the 
mean energy E m of the wave-packet ©, the greater is (At). For the present situation TPwkb is defined as [l9j . 



WKB 



2e J °- 



2^3(V-E)da 



/ 2 N /3(V - E)da 



(11) 



As an example, in Table HI we show a comparison between TP and TPwkb for different values of E m and Ai for 
the case with A = 0.001 and B = 0.001. We can see, clearly, that both tunneling probabilities coincide if we consider 
the appropriate At, for each E m . 



E m 


TP 


TPwkb 


At 


79 


5.74738 x 1CT 303 


8.88570 x 10~ 303 


12.5 


99 


1.22758 x 1CT 256 


3.03907 x 10~ 257 


18.5 


119 


8.00611 x 10~ 212 


2.78301 x 10~ 213 


26 


139 


9.09170 x 10~ 169 


1.31219 x 10~ 170 


33 


159 


7.49408 x 10~ 129 


5.04561 x 10~ 129 


45 


179 


1.78402 x 10~ 89 


2.22444 x 10~ ss 


57 


199 


3.42155 x 10~ 48 


1.46158 x 10~ 4S 


73 


209 


6.57786 x 10~ 29 


6.20233 x 10~ 29 


82 


216 


8.14639 x 10~ 16 


2.67764 x 10" 15 


89 


222 


5.65591 x 10 -03 


1.11987 x 10~ 03 


100 



TABLE I: A comparison between TP and TPwkb for 10 different values of E m with its associated integration time Ai for the case with 
A = 0.001 and B = 0.001. 
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FIG. 3: Tunneling probability as a function of the mean energy of the initial wave packet for A = 0.001, B = 0.001. 

In order to have an idea on how TP may differ from the TPwkb, we let the initial wave-packet (|T0| . with different 
mean energies, evolve during the same time interval At. We consider the example given above, in this section, with 
a common time interval of 100. We show this comparison in Table [TTl in the appendix, where we have an entry for 
TPwkb- It means that, we computed the values of the TPwkbs for each energy used to compute the TPs, in the 
case where A = 0.001 and B = 0.001. As we can see from Table ITT1 for this choice of At the tunneling probabilities 
disagree for most values of E m . They only agree for values of E m very close to the top of the potential barrier. There, 
because the values of E m are similar to 222, At is almost sufficient to guarantee that \& r of each wave-packet does 
not contribute to the TP. 



B. Tunneling probability as a function of A and B 

The Chaplygin gas has two parameters A and B equation |(5J). In the present subsection, we study the behavior of 
TP equation (TIT)]) as a function of both parameters. It is important to notice from the Chaplygin potential, equation 
([2]), that the parameter A behaves as a increases very much like a positive cosmological constant. Therefore, we 
expect that TP grows with A in the same way as it grows with the cosmological constant [ll[. On the other hand, 
the dependence of TP on the parameter B can not be inferred from previous works. 

We start studying the dependence of TP on the parameter A. In order to do that, we must fix the initial energy E m 
for the radiation and the parameter B. Then, we must compute the TP for various values of the parameter A. From 
that numerical study we conclude that the tunneling probability grows with A for fixed E m and B. As an example, 
we consider 21 values of A for fixed E m = 187 and B = 0.001, such that, the maximum energy of the potential barrier 
(PE max ), for each A, is greater than 187. The values of A, TP, a\ and are given in table HTT1 in the appendix. 
With those values, we construct the curve TP versus A, shown in figure [4] Due to the small values of some TP's, 
we plot the logarithms of the TP's against A. From figured! it is clear that the tunneling probability increases with 
A for a fixed E m and B, as expected. Therefore, it is more likely for the universe, described by the present model, 
to nucleate with the highest possible value of the parameter A. That result is in agreement with the approximated, 
analytical expression for TP as a function of A, obtained in Reference (T3 |. with the boundary conditions of the 
tunneling wave-function [l3| . 

Let us consider, now, the dependence of TP on the parameter B. In order to do that, we must fix the initial energy 
E m for the radiation and the parameter A. Then, we must compute the TP for various values of the parameter 
B. From that numerical study we conclude that the tunneling probability grows with B for fixed E m and A. As 
an example, we consider 21 values of B for fixed E m = 222 and B = 0.001, such that, the maximum energy of the 
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FIG. 4: Tunneling probability as a function of the parameter A for E m = 187 and B = 0.001. 



potential barrier (PE max ), for each B, is greater than 222. The values of B, TP, a% and ai are given in table iTVl in 
the appendix. With those values, we construct the curve TP versus B, shown in figure [3 Due to the small values 
of some TP's, we plot the logarithms of the TP's against B. From figure El it is clear that the tunneling probability 
increases with B for a fixed E m and A. Therefore, it is more likely for the universe, described by the present model, 
to nucleate with the highest possible value of the parameter B. That result is in agreement with the approximated, 
analytical expression for TP as a function of B, obtained in Reference [l2^ . with the boundary conditions of the 
tunneling wave-function [l3[ ■ 
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FIG. 5: Tunneling probability as a function of the parameter B for E m = 222 and A = 0.001. 
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VI. INITIAL CONDITIONS FOR INFLATION 



Based on these results one can evaluate the duration of the Planck era according to these models. In other words, 
the time that the universe takes to emerge at the right side of the potential barrier. According to the literature [20| . 
the time that a particle takes to cross a potential barrier is given by 



planck 



2a x 
TP' 



(12) 



where a\ corresponds to the turning point to the left of the barrier and TP is the tunneling probability Eq. f) 1 0[) . We 
have computed t p i anc k, for all examples considered here. In Tables HT1 IIIII and IIV1 in the appendix, we have an entry 
to t p i 

After the time interval t p i anc k the universe emerges to the right of the barrier with a 2 and p a = 0, because a 2 
is a turning point. Now, one may use these values as initial conditions to compute the dynamical evolution of the 
scalar factor through the integration of the appropriate classical Hamilton's equations. Due to our potential the 
classical evolution of a will correspond to an inflationary phase of the model, posterior to the Planck era. In the 
inflationary phase, according to the literature [211 ] . the universe has increased its size in many orders of magnitude, 
in a very short time. The table HH in the appendix, shows the initial conditions (a 2 ) for different values of E m , where 
the corresponding moments p a always vanish. Applying the 8th order Runge-Kuta method, implemented through 
a numerical routine in MAPLE, we have integrated the Hamilton's equations for many sets of initial conditions, in 
order to calculate the time evolution of the scale factor (conformal time) . We have concluded from that numerical 
investigation that the scale factor always increases many orders of magnitude in very short period of time, which 
confirms the development of a inflationary phase just after the Planck era. 

As an example, consider the case where the wave packet has E m — 220 and initial conditions a 2 = 12.95147904 
and p a = at r = 0, the scale factor increases from 12.95147904 to 11465.3370549124 in the interval At = 1.46. 
This behavior is shown in figure [6l which presents (in logarithm scale) the time evolution of the scale factor a(r), as 
a function of the conformal time. 



FIG. 6: Scale factor as a function of the conformal time (logarithmic scale). The initial conditions are 02 = 12.95147904 and 
p a = in r = 0, after the Planck era for wave packets with E m = 220. 



VII. CONCLUSION AND FINAL REMARKS 



The quantization of a FRW model with radiation and Chaplygin gas was made via finite differences method and 
the Crank-Nicholson scheme. We considered initial wave packets with finite norm, and well defined values of the mean 
energy, associated to the radiation fluid, and the parameters of the Chaplygin gas. Using this numerical method, one 
builds the time evolution of these wave packets. The results show that the quantum dynamics of these models may 
be described by wave packets with finite norm, which satisfy appropriate boundary conditions. In other words, such 
wave packets are well-defined over the whole space, even when the scale factor goes to zero. 
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Observing these wave packets, one can see the presence of strong oscillations for the values of the scale factor, 
which precede the effective potential barrier of the model. Then, when the packets reach the potential barrier, the 
amplitude of those oscillations decrease drastically, vanishing at a point far away from the barrier, at infinity. This is 
the quantum tunneling, which explains the birth of the universe with a certain, well-defined, size after the tunneling. 
Via numerical calculations, it was possible to determine the dependence of the tunneling probability with E m , A 
and B. We also compared it with the tunneling probability obtained using the WKB approximation. The results 
showed that both increase with the mean energy of the initial wave packets and show a good agreement, under certain 
circumstances. We verified that the duration of the quantum phase depends on the mean energy assigned to the 
initial wave packet, i.e., the bigger the energy, the shorter is the duration of the Planck era. We also, determined the 
conditions of the universe at the end of the Planck era. These conditions, obtained by the quantization of the model, 
were used as initial conditions to determine the classical evolution of the scale factor of the universe. The evolution 
shows that, after its quantum phase, the universe undergoes a very accelerated growth. It increases its size in many 
orders of magnitude in a very short time interval. We identify this phase as the inflationary phase. 
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APPENDIX A: TABLES 



Energy (E m ) 


PT 


tplanck 


ai 


ai 


PTwKB 


222.000000 


0.00565591 


4134.35 


11.691761 


12.701983 


1.1198667720 • 10~ 3 


221.000000 


0.000963791 


23947 


11.539957 


12.840057 


1.307940728 ■ 10~ 5 


220.000000 


0.0001864 


122476 


11.414765 


12.951479 


1.521014497 ■ 10~ 7 


219.000000 


3.99923 • 10~ 5 


565372 


11.305263 


13.047172 


1.762093474 ■ 10~ 9 


218.000000 


9.35514 • 10~ 6 


2.39578 • 10 6 


11.206415 


13.132172 


2.033617858 ■ 10" 11 


217.000000 


2.35555 • 10~ 6 


9.43763 • 10 6 


11.115403 


13.209295 


2.338009646 ■ 10~ 13 


216.000000 


6.32225 • 10~ 7 


3.48941 • 10 7 


11.030465 


13.280306 


2.677635300 ■ 10~ 15 


215.000000 


1.79509 • 10~ 7 


1.22004 • 10 8 


10.950407 


13.346395 


3.054754108 ■ 10~ 17 


214.000000 


5.3592 ■ 10~ 8 


4.05821 • 10 8 


10.874382 


13.408411 


3.471460808 ■ 10~ 19 


213.000000 


1.674 ■ 10~ 8 


1.29054 • 10 9 


10.801756 


13.466987 


3.929622371 ■ 10~ 21 


212.000000 


5.4483 ■ 10~ 9 


3.9396 ■ 10 9 


10.732048 


13.522604 


4.430810920 ■ 10~ 23 


211.000000 


1.84123 • 10~ 9 


1.15845 • 10 10 


10.664877 


13.575643 


4.976232532 ■ 10~ 25 


210.000000 


6.44182 • 10~ 10 


3.29098 • 10 10 


10.599935 


13.626410 


5.566654752 ■ 10~ 27 


209.000000 


2.32726 • 10~ 10 


9.05525 • 10 10 


10.536974 


13.675154 


6.202333780 ■ 10~ 29 


208.000000 


8.66257 • 10~ 14 


2.41863 • 10 11 


10.475786 


13.722084 


6.882942312 ■ 10~ 31 


207.000000 


3.31557 • 10~ 14 


6.28321 • 10 11 


10.416196 


13.767372 


7.607503240 ■ 10~ 33 


206.000000 


1.30263 • 10~ 14 


1.59033 • 10 12 


10.358055 


13.811169 


8.374327820 ■ 10~ 35 


205.000000 


5.24514 • 10~ 12 


3.92792 • 10 12 


10.301236 


13.853600 


9.18096264 • 10~ 37 


204.000000 


2.16153 • 10~ 12 


9.47999 • 10 12 


10.245627 


13.894776 


1.002414945 ■ 10~ 38 


203.000000 


9.10506 • 10~ 13 


2.23856 • 10 13 


10.191134 


13.934793 


1.089979620 ■ 10~ 40 


202.000000 


3.91585 • 10~ 13 


5.17776 • 10 13 


10.137670 


13.973737 


1.180296533 ■ 10~ 42 


201.000000 


1.71766 • 10~ 13 


1.17429 • 10 14 


10.085160 


14.011682 


1.272788104 ■ 10 -44 


200.000000 


7.67719 • 10~ 14 


2.61386 • 10 14 


10.033537 


14.048694 


2.059920186 ■ 10~ 46 


199.000000 


3.49333 • 10~ 14 


5.71532 • 10 14 


9.982742 


14.084834 


1.461582575 ■ 10~ 48 


189.000000 


3.05503 • 1CT 17 


6.22616 • 10 17 


9.510554 


14.407881 


2.255557889 ■ 10~ 68 


179.000000 


8.61288 • 10~ 20 


2.10926 • 10 20 


9.083391 


14.680929 


2.224440858 ■ 10~ 88 


169.000000 


5.67252 • 10~ 22 


3.06189 • 10 22 


8.684327 


14.920460 


1.364184268 ■ 10~ 108 


159.000000 


7.12611 • 10~ 24 


2.33052 • 10 24 


8.303768 


15.135558 


5.045605144 ■ 10~ 129 


149.000000 


1.48584 • 10~ 25 


1.06815 • 10 26 


7.935486 


15.331854 


1.087228346 ■ 10~ 149 


139.000000 


4.64724 • 10~ 27 


3.26 • 10 27 


7.575013 


15.513119 


1.312194871 ■ 10~ 170 


129.000000 


2.01877 • 10~ 2S 


7.15172 • 10 28 


7.218847 


15.682026 


8.478448344 ■ 10~ 192 


119.000000 


1.14576 • 10~ 29 


1.19816 • 10 30 


6.864019 


15.840547 


2.783009253 ■ 10~ 213 


109.000000 


8.07577 • 10~ 31 


1.61169 • 10 31 


6.507811 


15.990187 


4.364145700 ■ 10~ 235 


99.000000 


6.76014 • 10~ 32 


1.81876 • 10 32 


6.147539 


16.132125 


3.039077147 ■ 10~ 257 


89.000000 


6.44132 • 10~ 33 


1.79477 • 10 33 


5.780359 


16.267303 


8.604523368 ■ 10~ 280 


79.000000 


6.68111 • 10~ 34 


1.61741 • 10 34 


5.403039 


16.396490 


8.885695860 ■ 10 -303 


69.000000 


7.14821 • 10~ 35 


1.40222 • 10 35 


5.011661 


16.520320 


2.918766558 ■ 10~ 326 


59.000000 


7.31694 • 10~ 36 


1.25767 • 10 36 


4.601147 


16.639326 


2.553230826 ■ 10~ 3SO 


49.000000 


6.33942 • 10~ 37 


1.31383 • 10 37 


4.164443 


16.753960 


4.679948348 ■ 10~ 375 


39.000000 


3.66812 • 10~ 3S 


2.01243 • 10 38 


3.690908 


16.864607 


1.278053730 ■ 10- 400 


29.000000 


8.13428 • 10~ 40 


7.77617 • 10 39 


3.162678 


16.971602 


1.801636092 ■ 10 215 


19.000000 


1.21515 • 10~ 42 


4.18787 • 10 42 


2.544451 


17.075236 


2.645204325 ■ 10~ 455 


9.000000 


5.24188 • 10~ 51 


6.64299 • 10 50 


1.741086 


17.175763 


1.109066979 ■ 10~ 485 


8.000000 


1.28841 • 10~ 51 


2.54671 • 10 51 


1.640607 


17.185653 


6.395403152 ■ 10~ 489 


7.000000 


8.60797 • 10~ 52 


3.56372 • 10 51 


1.533820 


17.195515 


3.262847086 ■ 10~ 492 


6.000000 


4.88772 • 10~ 52 


5.80764 • 10 51 


1.419304 


17.205349 


1.446721601 ■ 10~ 495 


5.000000 


3.13507 • 10~ 52 


8.26144 • 10 51 


1.295012 


17.215155 


5.443204796 ■ 10~ 499 


4.000000 


1.78021 • 10~ 52 


1.30074 • 10 52 


1.157802 


17.224932 


1.680368783 ■ 10~ 502 


3.000000 


9.70346 • 10~ 53 


2.06603 • 10 52 


1.002384 


17.234682 


4.045122584 ■ 10~ 506 


2.000000 


4.42249 • 10~ 53 


3.70114 • 10 52 


0.818412 


17.244404 


6.963887216 ■ 10~ 510 


1.000000 


1.31081 • 10~ 53 


8.8352 • 10+ 52 


0.579062 


17.254099 


7.140705628 ■ 10~ 514 



TABLE II: Relation between the mean energies of the wave packets and the tunneling probability for A = 0.001; B = 0.001; -/V = 4000 
(spacial discretization); a max = 30; tmax = 100; and V m ax = 223.5282023, which corresponds to the maximum value of the potential. 
Here, we see as well the initial conditions a2, for the Universe to emerge from the Planck era to the Inflationary phase, all of them with 
Pa = 0. 
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A 


TP 


tplanck 


Oi 


«2 


0.001000 


8.70789o-18 


2.16406o+18 


9.422222 


14.465801 


0.001020 


2.77614o-15 


6.81336o+15 


9.457412 


14.340803 


0.001040 


5.7598o-13 


3.2965o+13 


9.493610 


14.216939 


0.001060 


7.76552o-ll 


2.45467c+ll 


9.530909 


14.094026 


0.001080 


6.76155o-09 


2.83054o+09 


9.569412 


13.971873 


0.001100 


3.75605c-07 


5.11668c+07 


9.609242 


13.850280 


0.001120 


1.30432c-05 


1.47978c+06 


9.650536 


13.729031 


0.001140 


0.00027424 


70693.3 


9.693459 


13.607891 


0.001160 


0.00332339 


5860.41 


9.738201 


13.486603 


0.001180 


0.0215238 


909.224 


9.784991 


13.364874 


0.001200 


0.0668755 


294.102 


9.834106 


13.242367 


0.001220 


0.0946396 


208.916 


9.885885 


13.118686 


0.001240 


0.105477 


188.492 


9.940751 


12.993352 


0.001260 


0.125288 


159.62 


9.999247 


12.865774 


0.001280 


0.139881 


143.866 


10.062084 


12.735189 


0.001300 


0.149677 


135.361 


10.130225 


12.600590 


0.001320 


0.162884 


125.305 


10.205035 


12.460569 


0.001340 


0.175986 


116.925 


10.288549 


12.313046 


0.001360 


0.187801 


110.586 


10.384061 


12.154690 


0.001380 


0.19885 


105.583 


10.497560 


11.979472 


0.001400 


0.212192 


100.309 


10.642415 


11.773989 



TABLE III: The computed values of TP, A, < p ; anc fc, a\ and a,2 for 21 different values of A when E = 187 and B = 0.001. 



B 


TP 




ai 


0,2 


0.001000 


0.00565591 


4134.35 


11.691761 


12.701983 


0.001350 


0.00565618 


4134.16 


11.691764 


12.701982 


0.001700 


0.00565644 


4133.97 


11.691766 


12.701980 


0.002050 


0.0056567 


4133.78 


11.691768 


12.701978 


0.002400 


0.00565696 


4133.59 


11.691770 


12.701977 


0.002750 


0.00565721 


4133.4 


11.691773 


12.701975 


0.003100 


0.00565746 


4133.22 


11.691775 


12.701973 


0.003450 


0.00565771 


4133.04 


11.691777 


12.701971 


0.003800 


0.00565796 


4132.86 


11.691779 


12.701970 


0.004150 


0.00565821 


4132.68 


11.691781 


12.701968 


0.004500 


0.00565846 


4132.5 


11.691784 


12.701966 


0.004850 


0.00565871 


4132.32 


11.691786 


12.701964 


0.005200 


0.00565895 


4132.14 


11.691788 


12.701963 


0.005550 


0.00565919 


4131.96 


11.691790 


12.701961 


0.005900 


0.00565944 


4131.79 


11.691793 


12.701959 


0.006250 


0.00565968 


4131.61 


11.691795 


12.701958 


0.006600 


0.00565992 


4131.43 


11.691797 


12.701956 


0.006950 


0.00566016 


4131.26 


11.691799 


12.701954 


0.007300 


0.0056604 


4131.08 


11.691801 


12.701952 


0.007650 


0.00566064 


4130.91 


11.691804 


12.701951 


0.008000 


0.00566088 


4130.74 


11.691806 


12.701949 



TABLE IV: The computed values of TP, B, t p i anc k, a\ and <J2 for 21 different values of B when E = 222 and A = 0.001. 



